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Abstract 

We show that one can obtain asymptotic evolving boost-invariant 
geometries in a simple manner from the corresponding static solutions. 
We exhibit the procedure in the case of a supergravity dual of R- 
charged hydrodynamics by turning on a supergravity gauge field and 
analyze the relevant thermodynamics. Finally we consider turning on 
the dilaton and show that electric and magnetic modes in the plasma 
equilibrate before reaching asymptotic proper times. 



1 Introduction 

Quark-gluon plasma, the deconfined state of matter produced at RHIC is 
currently the focus of numerous experimental and theoretical studies. Since 
it appears that the gauge theory coupling in the plasma is large P, it is 
interesting to develop nonperturbative methods for studying its properties. 

A very effective framework for studying nonperturbative phenomena is 
the AdS/CFT correspondence |2j|. Even in its simplest form, for Af = 4 
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Super- Yang-Mills theory, at finite temperature one may expect many simi- 
larities with properties of strongly coupled deconfined QGP. 

A lot of work has been done in the study of properties of static plasma 
at fixed temperature e.g. transport coefficients have been computed [3lll], 
the drag force acting on a moving quark was investigated [5j. Much less 
is known regarding questions related to time- dependent phenomena like the 
origin of thermalization, hydrodynamic evolution etc. Works investigating 
these issues include [HI [3 El El [lOl [11] . 

In order to gain more understanding of the above mentioned dynamical 
processes in strongly interacting gauge theory using AdS/CFT it is worth- 
while to investigate the structure of dual time-dependent geometries. In 
general holographic renormalization [12] gives a prescription for construct- 
ing a dual geometry to any given gauge theory energy-momentum spacetime 
profile. In ^ it was advocated that the requirement of nonsingularity of the 
dual geometry selects the physical energy-momentum evolution in the gauge 
theory. At leading order, in a boost invariant setting, this requirement led to 
late-time perfect fluid hydrodynamical evolution [S], while carrying out the 
analysis also for subleading asymptotic time^ [TT] determined the effects of 
viscosity on the evolution with the exact (shear) viscosity coefficient. 

The dual geometry for late asymptotic times constructed in [8] bears a 
remarkable similarity to the static black hole but with the position of the 
horizon moving with a specific scaling with proper-time. In this note we 
would like to perform an analogous construction for the more complicated 
case with R-charged matter, show that a similar phenomenon also occurs in 
this case and exhibit the origin of such a behavior. 

Another direction of generalization of the evolving geometry is to turn on 
the dilaton field. This has the physical interpretation on the gauge theory 
side of allowing for differing expectation values of squares of electric and 
magnetic fields in the evolving plasma. We find that, for asymptotic proper- 
times, in order to have a nonsingular geometry one has to have equilibration 
between electric and magnetic modes. 

The plan of this paper is as follows. First in section 2 we will briefly 
review boost-invariant hydrodynamics with a conserved charge. In section 3, 
we will construct the dual charged evolving geometry for asymptotic proper- 
times and explain the origin of its marked similarity to the corresponding 
static geometry. We will then examine, in section 4, the thermodynamics of 

^And using some results of [TO] , 
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the resulting evolving system. Finally, in section 5, we will consider turning 
on the dilaton. We close the paper with concluding remarks. 



2 Hydrodynamics with a conserved charge 

Perfect fluid hydrodynamics with a conserved charge is described by the 
energy-momentum tensor and the current 

T''-' ={e + p)uV + pri'"' , = pu^" , (2.1) 

which are conserved i.e. 

a^T^^ = 0, d^r = 0. (2.2) 

Moreover for the conformal theory that we are considering T'*^ = which 
gives e = 3p. 

Let us now impose the requirement of boost-invariance in the longitudinal 
plane and no dependence on transverse coordinates. In the natural propcr- 
time/spacetime rapidity coordinates the Minkowski metric has the form 

ds^ = -dr^ + T^dy^ + dx^ . (2.3) 

Boost-invariance then forces the fluid velocity to be — (1,0,0,0), and 
energy momentum conservation leads to Bjorken evolution, 

<r)^^, (2.4) 

T3 

while current conservation leads to 

V^J^^^ + -r = J'' = (^, 0,0,0). (2.5) 

dT T T 

We will show below that this scaling leads to a nonsingular dual geometry. 

3 Supergravity analysis 

The 5D Einstein Maxwell action for the the minimal cases is given by 

I = j [y^{R + 12 - \FapF^^) - ae^^^'^'F^f^F^^A,^ . (3.1) 
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where G5 = tx/{2NI) and a = 1/(48^/3) in our convention. Any solution 
of the above action can be consistently embedded into the lOD type IIB 
supergravity [I3j. The gravity equations of motion become 

1^1 

Ra/3 = —4:gal3 — J^F^uF^" Qap + -Fa ^Fp^ , (3.2) 

VaF"'^ - -^e^'^^^^FapF^s = . (3.3) 

The static black hole solution dual to the finite temperature A/" = 4 SYM 
plasma with R-charge chemical potential is given by 

ds^ = \(-h(z)dt^ + dx^ + dzyh(z)) , (3.4) 
z^ 

Fzt = qz , (3.5) 

where h{x) = 1 — ax^ + and all other components of the gauge field 

are vanishing. 

In order to analyze the gauge theory energy-momentum tensor in the 
most convenient manner one usually passes to the Fefferman-Graham coor- 
dinates [12]. However it turns out that contrary to the uncharged case g = 
considered in [8], the Fefferman-Graham form of (13. 4p cannot be expressed 
in terms of elementary functions. Therefore it is more convenient to adopt 
coordinates similar to (13.41) also for the evolving case. 

Namely the proper-time dependent metric ansatz will be 



F,. = K{z,t) , (3.7) 



which is a general ansatz compatible with boost invariance. Note that with 
respect to the Fefferman-Graham form of coordinates, the z and r coordinates 
have to be redefined. With this ansatz, the gauge field equation part can be 
solved with an integration constant q by 



Z A+D~B 

Fzr = q-e^- . (3.8) 



Then the scalar F^ becomes 



^6 

= -2q^—e-^ . (3.9) 
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To study nontrivial late time scaling behaviors, we take the scaling variable 
as f = z/t^/'^ and work in the large r limit while keeping v fixed. Then in 
the scaling limit, one has 

A = a{v) +0{1/t^) , B=b{v) +0{1/t'^) , D=d{v) + 0{1/t^) , {'i.l^)) 
F' = f\v) + 0{l/T^) , (3.11) 

where \] denotes some positive power so that the terms can be ignored in 
the large proper-time limit. One can, of course, consider the case where 
= but then the analysis will be reduced to the problem of uncharged 
fiuid dynamics. 

From (13.91) . one concludes then that, in the scaling variable v = the 
power s should be fixed as 4/3 once we have a nonvanishing charge density 
in the scaling limit. 

Let us rewrite the Einstein equation (13. 2p in the following form 

R'^P = (4 + ^F'^ + ^F-^Fp, . (3.12) 

The scaling large r limit of the diagonal equations (rr, zz, yy and xx com- 
ponents) is 

R^r = -4 + R', = -4 + ^F^ (3.13) 

R\ = -4 - ^F^ R\ = -4 - , (3.14) 

where the explicit formulas for the R'^p in the scaling limit are given in the 
appendix. Apart from these equations we also have the leading part of the 
off-diagonal equation K^^ = 0: 

v{2b"{v)+b'{vf)+2{3a'{v)+3d'{v)-2b'{v))-v{a'{v)+d'{v))b'{v) = 0. (3.15) 

Remarkably enough a direct analog of the static solution: 



a{v) = \ogil-av^ + ^v^] (3.16) 



t 
12 

b{v) = (3.17) 
d{v) = -log(l-av^ + ^vA (3.18) 



solves all the above equations! If we recall now the definition of the scaling 
variable v = z/t^^^, we see that asymptotically for large r the geometry looks 
like the static charged black hole with 

h = l-a^ffir)z' + ^^il^z' (3.19) 

but with the parameters a^jfij) and qeffij) being r-dependent. In particular 
the scaling is exactly such that the effective charge behaves like 

qeffir) ~ ^ (3.20) 

as follows from (12.51) . We shall return to the more precise determination of 
charge density later on. 

Before we examine in detail the physical properties of the evolving solu- 
tion let us try to understand why such a simple generalization of the static 
solution gives an evolving boost invariant solution. 

It turns out that the diagonal components of the Ricci tensor i?"^ com- 
puted for the static metric 

ds^ = 4 i-e^^'^dt^ + e^^'^dy"^ + e<'^dx\ + e'^^'^dz'^) (3.21) 
z^ 

coincide with the scaling limit of diagonal components of the Ricci tensor 
R^jj of the boost-invariant evolving metric 

z^ 

with the substitution of z for v = z/t^^^. A similar scaling property holds 
here for the gauge field, hence equations (I3.13I) -( 13TT^ for the evolving case 
coincide with the corresponding equations for the static metric with v inter- 
changed with z. 

Once equations (I3.13p - (l3.14p are solved, the off-diagonal equation (13.151) 
is automatically satisfied without giving any further independent restriction. 
This can be shown as follows. Note the Bianchi identity, 

V^i/n. = , (3.23) 
where H'^^ is the tensor appearing in the usual form of the equation of motion 

H^', = - ]^6^R - S-kG^T^ = (3.24) 
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with including the contributions from the gauge fields and the cosmo- 
logical constant. The equation fl3.23p is a mathematical identity once the 
solution fl3.8p of the Maxwell equations is used, which insures the covari- 
ant conservation of the tensor T^. Then v = t component of the Bianchi 
identity can be written as 

d,{V^H%) + dr{^/^H\) - v^r^^if^, = . (3.25) 

2 

Note that H\ \^ 0{z/t) {= 0{v /t^)) in its leading order. Assuming (I3.13P - 

2 

fl3.14p are solved in their leading order, if^^ = H^y = H\ = 0{v/ts) at most. 
The terms involving H^^ = H^y = H\ in fl3.23p have an extra r derivative 

2 

and, therefore, are w/rs higher order than the leading terms of H'^t-. Hence 
by setting those higher order terms to zero, one is left with 

d,{y/^H\)=Q. (3.26) 

Fixing the integration constant to zero by the asymptotic condition, one does 
have the relation H\ = R\ = 0. Thus the off-diagonal equation automati- 
cally follows from the remaining part of the equations. 

Let us emphasize that the above procedure gives a way to obtain the form 
of asymptotic geometry only for large proper-times from the corresponding 
static solution. The resulting solution is, however, not an exact solution for 
smaller r. The sub-asymptotic form of the metric does not seem to be any 
longer linked to the static solution. In a way this is not surprising since the 
nonlinear evolution includes effects of viscosity and even more pronounced 
deviations from perfect fluid hydrodynamics for small proper-times. What 
seems to be more surprising is the fact that any such correspondence between 
static and evolving boost-invariant solutions exists at all. 

4 Thermodynamics 

Let us now analyze the thermodynamic properties of the evolving geometry 
derived in the previous section. 

Let us flrst evaluate the entropy per unit volume. The boundary volume 
element here is given by 

dV = TdydxUx^ . (4.1) 

We will use just the extrapolated static formulas since in any case the asymp- 
totic geometry is not valid at the horizon but rather in the scaling limit 
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r — s> oo with V fixed. The location of the extrapolated horizon is thus ob- 
tained by finding a smaller positive root of 

l-avi + ^vl = Q, (4.2) 

so the entropy density is given by 

S N'^ I N"^ I 
^ = T7 = ^--.=^-7JZ (4-3) 



while the temperature is 



V 2tt z'l 2tt Vyj 



T = ' ^ = ^ r~3 . (4.4) 

In order to find the energy density e(r) one has to pass to Fefferman- 
Graham coordinates defined by 



and read off eir) from 



^(^) = TrT l™n 4 • (4-6) 



We have to redefine the z variable as 



z = z,g[1-i'-^ + ...] (4.7) 



and requiring 

-1 



z 



2 



(4.8) 

^"FG 

fixes the coefficient ^ = aj%. At this order and in the scaling limit one does 
not need to redefine r. Applying the redefinition (14. 7p to the g^-T component 
of the metric we use (14. 6 p to get 

e(T) = — ^ ^ = — ^ ■ -a ■ r s = ^— r s . (4.9) 

2vr2 rl 2vr2 4 8vrM(l - ^^D 



By the similar procedure for the other component of the metric, one finds 
that p = e/3. 

Finally the charge and chemical potential can be found from the behavior 
of the gauge field. Let us first work out the charge density. The displacement 
is evaluated as 

which is constant. The boundary charge density is then defined by 

p= dydx^dx^D'^V = -^q/r (4.11) 

showing the expected r dependenc^. The chemical potential can be given as 
the difference between the horizon and the boundary values of the Coulomb 
potential 

^^ = A,{zn)-A,{z = Q) = ]-■^. (4.12) 
One may check that the Gibbs potential Vl satisfies the relation, 

Vl/V = -p = e-Ts- lip. (4.13) 
Also the energy density can be written as 

e p) = — — — T 1 + ^r^) , 4.14 

from which one may check that the chemical potential fi is indeed conju- 
gated to p by the relation yU = de/dp. Finally there is the requirement of 
thermodynamic stability leading to the condition [TU [151 [IB] , 



^The boundary charge density and the chemical potential can be determined by study- 
ing the boundary behavior of the gauge fields using AdS / CFT dictionary. The results are 
the same as the ones from the bulk method used here. 
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5 The dilaton and electric/magnetic equili- 
bration 



It is interesting to consider turning on the dilaton instead of the gauge field. 
Such a setup would correspond to considering configurations (states) in gauge 
theory with a nonvanishing expectation value of tr F^^F'^'^ (note that this is 
now the 4D gauge theory field strength and not the 5D supergravity field 
considered before in this paper). This means that such a configuration of the 
plasma has 

'tiEA ^ (tiB^) . (5.1) 



The corresponding 5D action in the Einstein frame is given by 
1 r . /„ 1 



y v^(i? + 12 - ^g-^d^<j)dp(p) . (5.2) 



levrGs 

The supergravity equations of motion becomes 

Ra/3 = -4:gaf3 + ^dc,4>di3(f) , (5.3) 
V^0 = O. (5.4) 

As in Ref. [H], we shall use the Fefferman-Graham form of coordinate by 
setting D = this time while allowing a nontrivial C{z,t) = c{v) + 0{1/t^) 
with the general scaling variable v = z/r^^^ with < s < 4. Also in the 
boost invariant scaling limit, the scalar field behaves 

4>{z,t) = ^{v) + 0{1/?). (5.5) 

The scalar field equation in the leading order can be integrated leading 

to 

= ^W-; ■ (5-6) 
Using this expression, the remaining Einstein equations become 



-2v{a' + h' + 2c') + v\{a'f + {h'f + 2{c'f + 2(a" + h" + 2c")) 

= -2;^ ^ (5-7) 
-2t;(4a' + 6' + 2c') + t;'((a')' + a'(6' + 2c') + 2a") = , (5.8) 
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-2v{a' + Ah' + 2c') + v^{{h'Y + h'{a' + 2c') + 26") = , (5.9) 
-2v{a' + h' + 5c') + v^{2{c')^ + c'(a' + h') + 2c") = , (5.10) 

2t;(li^lz^ _ 6' _ 2c') - v\{h')^ - a'ih' + 2c') + 2idf + 26" + 4c") 
s 

= 2^. (5.11) 

Interestingly these equations can be solved exactly. The steps are as follows. 
First we add ( 15. 7p and ( 15. lip to cancel out the scalar contributions. We then 
linearly combine the resulting expression with (15. 8p such that all the second 
derivatives cancel out. The final result becomes 

(4 - 3s)a' + (s - 4)6' + 2sc' = . (5.12) 

Using the boundary condition a(0) = 6(0) = c(0) = 0, this is solved by 

a = M-2m, b = M+{2s-2)m, c = M + {2 - s)m . (5.13) 

with M(0) = m(0) = 0. Inserting these expressions into Eqs. fl5.8p - fl5.10p . 
one finds that the three equations are reduced to 

-7vM' + v^{2{M')^ + M") = 0, (5.14) 

-3vm' + v^{2m'M' + m") = . (5.15) 

These are solved by [S] 

M = ^ ln(l - A^t;^) , (5.16) 

m = i-lnf ^-f' ') . (5.17) 

What remains is to satisfy (15. 7p with the above expression of M and m, 
which leads to 

3.^_8. + 8^ (5,18) 

24 96 ^ ^ 

Then finally computing Rf^^^pRi""^^ [8], one can check that it is nonsin- 
gular only for k = and s = 4/3. Therefore one concludes that the evolving 
geometry with a nontrivial scalar field is not allowed for the late time. 

The above result shows that fiuctuations of electric and magnetic modes 
tend to equilibrate very fast. Perhaps faster than power-like behavior but 
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note that above result only shows that the scalar field should approach zero 
faster than with t] > for the finite scaling coordinate v. 

This is analogous to the stability property of the evolving perfect fiuid ge- 
ometry with respect to small scalar perturbations |9] in the linearized theory. 
These are in effect quasinormal modes which exhibit exponential decay. The 
above analysis extends this result to states far away from electric/magnetic 
equilibration. 

6 Conclusions 

In this note, we obtain the asymptotic evolving boost-invariant geometry in- 
volving conserved R-charge. Thereby the boost-invariant late time dynamics 
of strongly coupled A/" = 4 Super Yang-Mills with R-charge turned on are 
studied via the AdS/CFT correspondence. The result shows that the boost 
invariant late time state has necessarily to be in the perfect-fluid hydrody- 
namic regime even including the R-charge. 

The asymptotic large proper-time boost-invariant evolving geometry can 
be seen to arise from the corresponding static solution by substituting the 
scaling variable v for z in the coefficient functions. This property arises due to 
the specific form of the Ricci tensor in the scaling limit which closely mirrors 
the Ricci tensor for the analogous static solution. Let us note, however, that 
the asymptotic geometry is not an exact solution of the Einstein equations 
and sub-asymptotic corrections exist. These corrections are important as 
they encode specific time-dependent dynamical effects. 

We also discuss the 5D Einstein-scalar theory to show that the late time 
boost invariant geometry with a nontrivial scalar field is not allowed. In the 
gauge theory side this implies that a difference between electric and magnetic 
modes does not survive into the asymptotic proper-time regime. 
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Appendix. Expressions for the scaling limit of R^p 

Here we quote expressions for the diagonal components of the Ricci tensor 
used in the main text: 

R\ = ~e-''^''\lQ + v\2a"{v) + a\vf)+2v{h\v) + d'{v)-Aa'{v)) + 
W{h\v)+d\v))a\v)) (a.l) 

R\ = ~e~^^''\lQ + v\2a"{v) + a\vf)+v^{2h"{v) + h'{vf) + 

-2v{a\v) + h\v) - Ad\v)) - v\b'{v) + a'{v))d'{v)) (a.2) 

R% = ~e-''^''\l6 + v\2b"{v)+b'{vf)-2v{a'{v)-d'{v) + Ab'{v)) + 

-v^{d'{v)-a'{v))b'{v)) (a.3) 

i?^^ = -^{8 + vd'{v) -va'{v) -vb'{v)) . (a.4) 
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